Cotunneling through a magnetic single-molecule transistor based on N@C6o 
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We present an experimental and theoretical study of a magnetic single-molecule transistor based 
on N@C6o connected to gold electrodes. Particular attention is paid to the regime of intermediate 
molecule-lead coupling, where cotunneling effects manifest themselves in the Coulomb-blockade 
regime. The experimental results for the differential conductance as a function of bias, gate voltage, 
and external magnetic field are in agreement with our analysis of the tunneling rates and provide 
evidence of magnetic signatures in single-N@C6o devices arising from an antiferromagnetic exchange 
interaction between the Ceo spin and the nitrogen spin. 

PACS numbers: 



Single molecules carrying a magnetic moment, such 
as Mni2, Feg and N@C6o, have been studied inten- 
sively for their outstanding quantum properties. Some of 
these molecules are characterized by a large spin and ex- 
hibit macroscopic quantum tunneling [1, 2] and quantum- 
interference effects [3]. The ability to initialize and ma- 
nipulate the quantum state of a single magnetic molecule 
could open the road to new strategies for high-density 
information storage, quantum computing and molecu- 
lar spintronics [4, 5]. However, only few experiments 
have succeeded in fabricating single-molecule transistors 
(SMTs) [Fig. 1(a)] based on magnetic molecules [6-10]. 
Furthermore, only two experiments [7, 10] on single mag- 
netic molecules in three-terminal geometry have suc- 
ceeded to obtain the intermediate molecule-lead coupling 
regime. In this regime [11], cotunneling processes mani- 
fest themselves in the spectra of excited molecular states 
for a fixed charge state of the molecule, allowing an eas- 
ier and more precise characterization of magnetic prop- 
erties of single-molecule devices compared to sequential 
tunneling. For Mni2 clusters, steps in dl/dV have been 
observed [7] , but the dependence of the cotunneling steps 
on an external magnetic field was not investigated. Re- 
cently, a characteristic zero-field splitting and its evolu- 
tion in an external magnetic field has been attributed to 
the signature of anisotropy in a Fe4 single-molecule mag- 
net [10]. The aim of our work is then to study exchange- 
coupling-dependent cotunneling features in a single mag- 
netic molecule. 

Previous studies have shown that it is not straightfor- 
ward to identify the magnetic properties of a SMT since 
molecular magnetism is usually destroyed during the de- 
vice fabrication [7] due to strong interactions with the 
environment [12]. To conserve the magnetic properties, 
one option would be to incorporate a molecule whose 
magnetic moment is retained. In N@C6o, the spin of 
the encapsulated nitrogen atom, which sits in the cen- 




FIG. 1: (a) Colorized electron beam microscopy photography 
of an SMT consisting of a gold nanowire over an AI/AI2O3 
gate after electromigration. (b) Schematic picture of N@Ceo- 
(c) Spins states of neutral and charged N@C6o- 

ter of the fullerene molecule [13], is protected by the Ceo 
cage [Fig. 1(b)]. An exchange interaction between the 
C60 spin and the nitrogen spin permits an indirect mea- 
surement of the magnetic properties of the nitrogen [9]. 

Although experimental studies investigating the sta- 
bility of paramagnetism in N@Ceo on a substrate have 
not been performed so far, it is known from X-ray and 
ultraviolet photoelectron spectroscopy that Ceo changes 
its charge state when being adsorbed on noble metal sur- 
faces [14, 15]. This may lead to the formation of N@Cgg 
ions, when the number of electrons, nc 60 , added to the 
C60 cage is 1, or N@CgQ ions when nc 60 = 2. Since co- 
valently functionalized N@C6o is stable [16-18], various 
ions of N@C6o are believed to be stable [19]. 

The total spin of neutral N@C 60 is S = 3/2 [Fig. l.(c)]. 
However, the spin of the singly-charged ion N@CgQ can 
be S = 1 or S = 2, depending on whether the Ceo spin 
Sc 60 = 1/2 is antiparallel or parallel to the nitrogen spin 
[Fig. 1(c)]. N@Cgo has a fixed ground state spin S = 3/2 
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since the dianion cage is diamagnetic due to Jahn-Teller 
distortion. 

Following theoretical predictions [20], transport mea- 
surements on SMTs with N@C6o in the weak-coupling 
regime showed signatures of an antiferromagnetic ex- 
change interaction between the Ceo spin and the nitrogen 
spin [9] , in contradiction with the ferromagnetic exchange 
coupling obtained, for singly and triply ionized anions of 
N@C6o, using multi-configurational self-consistent field 
methods [21]. Therefore, we first present results that 
confirm the interpretation of previous experiment [9] in 
the weak-coupling regime. As the field of molecular elec- 
tronics in general and of molecular spintronics in par- 
ticular has been plagued by a lack of reproducibility, it 
is essential that we have been able to reproduce the re- 
sults for the prototypical single-molecule transistor based 
on N@C6o using independently fabricated devices. How- 
ever, since we succeeded in reaching the intermediate 
molecule- lead coupling regime, our results go beyond 
the sequential-tunneling regime. In the second part, we 
present results in the cotunneling regime. In agreement 
with our theoretical analysis, our experimental results 
provide evidence for an antiferromagnetic exchange in- 
teraction between the Ceo spin and the nitrogen spin. 

The SMT was fabricated by using the electromigra- 
tion technique [22]. The experiment was carried out in 
a dilution refrigerator at a base temperature of 35 mK. 
Our procedure and details of the measurement system 
are described in Ref. [23]. 

Figure 2(a) shows a color scale plot of the differential 
conductance dl/dV as a function of bias voltage V s d and 
gate voltage V g on the entire accessible gate range, per- 
formed during a first run of measurements. We empha- 
size that we observe only one charge degeneracy point (at 
gate voltage Vg ), where the energy of the nc 60 = 1 state 
aligns with energy of the nc 60 = 2 state. Figure 2(b) 
shows results for dl/dV in the vicinity of V® . Analyz- 
ing the slopes of the stability diagram, we obtain C g : 
Cr : Cl = 1 : 3.09 : 2.96 for the ratio of the capac- 
itances, where C gi Cr and Cl denote the capacitances 
of the gate, the source, and the drain with respect to 
the molecule. Thus the effective gate voltage is reduced 
by a factor of C g /(C g + C R + C L ) = 0.142. From this 
we infer that the energy spacing between the relevant 
charge degeneracy points is at least larger than 230 meV 
and conjecture that the measurement was performed for 
a single N@C6o in the junction [24]. 

When the energy level spacing between the molecular 
states is large compared to the thermal energy, the hy- 
bridization energies Tr and Tl can be estimated by an 
analysis of the current amplitudes at different gate volt- 
ages [25]. We find Y L ~ 1.65 /ieV and Y R ~ 16.5 fieV 
(Tr/Tl — 10). The assumption of asymmetric molecule- 
lead couplings explains the different peak heights in 
dl/dV at negative and positive bias [Fig. 2(b)] and is 
taken into account in the theoretical treatment discussed 
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FIG. 2: Color scale plot of the differential conductance dl/dV 
as a function of bias voltage V s d and gate voltage V s (a) on 
the entire accessible V g range and (b) in the vicinity of V^ D . 
(c) Numerical results obtained from Eqs. (4) and (5) for Sd — 
-3.1625 eV, U = 3.0 eV, J = -0.4 meV, C 9 : C R : C L = 1 : 
3.09 : 2.96, and r R /F L ~ 10. 

below [Fig. 2(c)]. 

We compare the bias- and gate-voltage dependences 
of the measured differential conductance to numerical 
results obtained from a rate-equation approach. The 
system is described by a Hamiltonian of the form H = 
#n@c 6 o + Pleads + H u where 

#N@C 60 = ^d^Ceo + ^ n C 60 (^C 60 - 1) 

- JS C60 -S N -^flB(55 fl0 +55) (1) 
describes the molecular degrees of freedom, 

Pleads = ^2 S 

^ a k ^ a & a kcr (2) 

a—L,R kcr 

describes the conduction electrons in the leads, and 

H t = X^« a «kcA +C4^kcr) (3) 

a—L,R kcr 

describes the tunneling of electrons between the Ceo cage 
and the leads. Here, d\ creates an electron with spin a 
and energy in the LUMO of the Ceo molecule, while 
a akcr creates an electron in lead a with spin cr, momen- 
tum k, and energy e a k- Tuning the gate voltage V g allows 
one to shift the orbital energy Sd- The number and spin 
operator of the Ceo electrons are nc 60 = d\d G and 
S c 60 = Y^acr* d\{cr aa > j2)d a >, respectively, with cr the vec- 
tor of the Pauli matrices. The exchange interaction with 
the nitrogen spin Sn is denoted by J and the external 
magnetic field applied along the z axis by 5, where g is 
the ^-factor and /ig is the Bohr magneton. 

The leading contribution to the tunneling rates (cor- 
responding to sequential tunneling) is of second order in 
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H t and can be obtained from Fermi's Golden Rule, 



r: 



(seq) 
n— m' ,a 



+ [l-f(e n -e n ,-fi a )}\C:, n \ 2 y (4) 

Here n and n' denote the initial and final many-particle 
state of N@C6o, and / is the Fermi function. The matrix 
element with C° n , = (n\d a \n f ) denotes the net overlap 
of the initial and final states of the molecule. The typ- 
ical sequential-tunneling rate is given by T a = 27r£ a 2z/ a , 
where u a denotes the lead density of states (ft = 1). To 
mimick the finite life-time broadening of the sequential- 
tunneling peaks, we used an effective temperature of 
250 mK in Eq. (4). 

The next-to-leading-order contribution to the tunnel- 
ing rates (including cotunneling) is of fourth order in 
H t and can be computed using a T-matrix formalism 
[26, 27]. We find 
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■ t n " -r e + ir) 
x/(e - Ha) [l - f(e n -e n >+e- /i a /)] . (5) 

The divergences arising from the energy denominators 
are regularized by replacing rj by a finite life-time broad- 
ening r = + Tr. In the Coulomb-blockade regime 
far from the sequential-tunneling threshold, the Fermi 
functions suppress the integrand in Eq. (5) outside of a 
narrow energy interval. Replacing e in the energy denom- 
inators by its mean on these intervals, we can perform 
the integral. Close to the sequential-tunneling thresh- 
old, where this approximation becomes invalid, the co- 
tunneling contribution is overwhelmed by the sequential 
tunneling in any case. 

The good qualitative agreement between Figs. 2(b) and 
2(c) suggests that the Hamiltonian [Eqs. (l)-(3)] contains 
the crucial ingredients to explain both sequential and co- 
tunneling features revealed in the experiment. We find 
that a calculation under the assumption of J ~ —0.4 meV 
gives the best agreement with the experiment. 

In order to verify that the charge degeneracy point 
involves the two states nc 60 = 1 and nc 60 = 2, we in- 
vestigated the dependence of dl/dV on an external mag- 
netic field B. These measurements were performed for 
the same sample but during a different run. Results of 
dl/dV for a fixed gate voltage greater than Vg D are shown 
in Fig. 3. The conductance peaks at positive and neg- 
ative bias voltage first move apart until the magnetic 
field assumes a critical value. Then they change slope 
approaching each other again. Here the best fit to exper- 
imental data is obtained for a calculation that assumes 
J ~ — 0.3meV. The other parameters are also changed 
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FIG. 3: Color scale plot of the differential conductance dl/dV 
as a function of bias voltage and magnetic field B. (a) 
shows experimental results while (b) shows numerical results 
in the regime where the ground state is nc 60 = 2, with 
£d = -3.1722 eV, U = 3.0 eV, J = -0.3 meV, C 9 : C R : 
C L = 1: 3.29 : 3.43, and F R /F L ~ 6.7. 



somewhat, see the caption of Fig. 3. The microscopic 
mechanism behind the change in the exchange interac- 
tion is not yet clear. The characteristic pattern revealed 
in Fig. 3 allows us to deduce the charge and spin of the 
relevant molecular states to the left and to the right of the 
degeneracy point, in agreement with the results reported 
by Grose et at. [9]. 

We now turn to the cotunneling regime where trans- 
port is dominated by the coherent tunneling of electrons 
through the source-molecule-drain structure [28, 29]. 
This regime is characterized by a fixed molecular charge 
state and by the occurrence of steps in dl/dV. The po- 
sitions of these steps are independent of the molecular 
charging energies and the gate voltage but do depend on 
the magnetic excitation energies. As a consequence, mea- 
surements in this regime can be more easily interpreted 
than in the sequential-tunneling regime, where two dif- 
ferent charge states are involved. 

For N@CgQ , the ground state is a degenerate spin 
5 = 1 multiplet while the excited state, having an en- 
ergy J compared to the ground state, is a spin S = 2 
multiplet at zero magnetic field. Turning on a magnetic 
field leads to a Zeeman splitting of the spin multiplets as 
indicated in Fig. 4(a). Selection rules for the cotunnel- 
ing processes require lAS^I = or 1. This condition is 
used to determine all possible transitions that we expect 
to observe as steps in dl/dV (Fig. 4(b)). Our measure- 
ment presented in Fig. 4(c) is in qualitative agreement 
with our theoretical analysis presented in Fig. 4(d), ob- 
tained with the same exchange coupling J = — 0.3meV, 
confirming in an another tunneling process, the antifer- 
romagnetic exchange coupling between the Ceo spin and 
the nitrogen spin. 

The discussion for the case V 9 > V® , where the ground 
state is a doubly-charged N@CgQ anion, proceeds in an 
analogous way and is included in the supplementary in- 
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FIG. 4: (a) Zeeman diagram of the multiplets £ = 1 and 
S = 2 for N@Cg^ (nc 60 = 1). (b) Schematic showing the 
expected cotunneling steps, (c) Experimental and (d) numer- 
ical results, with J = — 0.3meV, of dl/dV for a fixed gate 
voltage V g < Vg . Colored lines indicate the positions of the 
magnetic excitations illustrated in (a) and (b). 

formation. The most striking difference is the presence 
of a cotunneling step at very low bias, for zero magnetic 
field [Fig. 2(b)]. This feature might be due to magnetic 
anisotropy of the molecular spin [10], caused by the spin- 
orbit coupling of the nitrogen, which we expect to be 
enhanced by the hybridization with the gold leads. 

In summary, we have presented an experimental and 
theoretical analysis of an SMT based on N@C6o- Our 
results in the sequential-tunneling regime are in agree- 
ment with the findings previously reported [9], and are 
a step towards achieving the high degree of reproducibil- 
ity necessary for progress in the field of molecular spin- 
tronics. Furthermore, we have investigated the cotunnel- 
ing regime, where transport is dominated by the coher- 
ent transfer of electrons through the molecular junction. 
Since measurements in this regime can be more easily 
interpreted than in the sequential-tunneling regime, the 
characteristic two-dimensional plots of dl/dV as a func- 
tion of bias voltage and magnetic field serve a character- 
istic fingerprint of the magnetic state. In particular, co- 
tunneling data obtained as a function of magnetic field 
corroborate the antiferromagnetic exchange interaction 
between Ceo spin and the nitrogen spin. 
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SAMPLES PRESENTATION AND COMPARISON OF THE RESULTS 



Single-molecule transistors (SMT's) based on N@Ceo contacted by gold leads were pre- 
pared by blow drying a dilute toluene solution of the N@Ceo molecule onto a gold nano-wire 
realized on an AI/AI2O3 back gate. Before blow drying the solution, the nano- wires were 
cleaned with acetone, ethanol, isopropanol and oxygen plasma. We then performed our 
electromigration procedure at 4 K. We tested 274 junctions and present in this section 3 
samples having an addition energy greater than 200 meV (Fig.S.l) 
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Fig.S. 1: Color scale plot of the differential conductance dl/dV as a function of bias voltage V s d 
and gate voltage V g : a) first run of measurements, b) second run of measurements, c) and d) 
other N@C6o SMT's in the weak coupling regime. 
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In Fig.S.l(a) and 1(b) we present first the stability diagram of the sample presented in the 
Letter, obtained at zero magnetic field during the first and the second run of measurements, 
respectively. The sample has been kept at cryogenic temperature (T < 4 K) during the two 
runs. However, changes in our electronic setup (different I-V converters, inversion of the 
ground connection) have probably generated a small electrical discharge inducing modifica- 
tions of some measurement features, as exhibited by the modification of the capacitances 
and tunneling rates extracted from a fine analysis of the stability diagram. For the first run, 
we obtain C g : C R : C L = 1 : 3.09 : 2.96, and Y R /Y L ~ 10, while for the second run we 
get C g : C R : C L = 1 : 3.29 : 3.43, and T^/r L ~ 6.7. The SMT experienced then some 
unexpected changes inducing a different electrostatic environment leading to a modification 
of the anti-ferromagnetic exchange interaction J. The evolution of J will be addressed in 
the next section. 

The striking difference between the two runs is observed at the charge degeneracy point. 
In Fig. S. 1(b), corresponding to the second run, we measure an opening of this latter. As it is 
not present in the first run as seen in Fig.S.l(a), we can rule out any correlation between this 
observation and the cotunneling step at low bias voltage observed for nc 60 = 2 (V g > Vg D ). 
No further measurements were performed in order to explain this feature, but a similar 
behavior was observed in another SMT based on a single N@C6o, as presented in Fig. S. 1(d). 

MAGNETIC FIELD SPECTROSCOPY IN THE COTUNNELING REGIME FOR 

n c 60 1 

In this section, we compare the magnetic field spectroscopy performed in the cotunneling 
regime, where nc 60 = 1 (V g < Vg D ), for the two runs of measurements. As depicted in Figure 
4 of the Letter, when the magnetic field increases, the different multiplet states are Zeeman 
split. For a magnetic field B < B C) the ground sate is and taking into account the 

selection rules, we observe four different transitions, until the magnetic field reaches the 
critical magnetic field B c for which the ground state becomes |2, +2). Only two transitions 
are then allowed, obeying \AS Z \ = or 1, as observed in Fig. S. 2(a). Unfortunately, if 
the critical field B c was experimentally accessible for the second run corresponding to J = 
—0.3 meV (Fig. S. 2(a)), no spin transition could be observed for the first run (Fig. S. 2(b)) as 
the value of J = —0.4 meV was too high. Indeed, as J depends exponentially on the distance 
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Fig.S. 2: Experimental differential conductance as a function of magnetic field for a fixed gate 
voltage V g < Vg D (nc 60 = 1) for (a) J = —0.3 meV (second run) and (b) J = —0.4 meV (first 
run). 

between the nitrogen and the Ceo cage, a small deformation can modify J significantly. 
More precisely, J is given by an exchange integral, which depends strongly on the electronic 
wave function of the Ceo LUMO. This wave function is not only changed by geometric 
deformation but also by electrostatics. For example, image charges in the metal electrodes 
draw it outwards towards the electrodes and away from the nitrogen. A modification of the 
environment can thus also lead to a variation of J. 

MAGNETIC FIELD SPECTROSCOPY IN THE COTUNNELING REGIME FOR 

n c 60 = 2 

In the Letter, we raised the question of the unexpected cotunneling steps at low voltage 
bias, obtained at zero magnetic field in the nc 60 = 2 region where the total spin of N@CgQ 
is S = 3/2. In Eq. (1) of the Letter, we assume no anisotropy. In this case, performing a 
magnetic field spectroscopy, we should observe one cotunneling step increasing linearly with 
magnetic field (Fig.S. 3(b)) corresponding to the transition from the |3/2, +3/2) state to the 
|3/2,+l/2) state as shown in Fig.S. 3(a). In Fig.S. 3(c), we present the dl/dV measurement 
as a function of magnetic field. We clearly observe that the cotunneling steps evolve linearly 
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Fig.S. 3: (a) Zeeman diagram of the multiplet S = 3/2 for N@Cgg ( n C 60 = 2). (b) Schematic 
showing the expected cotunneling steps deduced from the Zeeman diagram, (c) and (d) dl/dV 
measurements as a function of magnetic field for a fixed gate voltage. 

in the high magnetic field limit. 

However, when we focus on the low magnetic field spectroscopy measurements 
(Fig.S. 3(d)), we observe a non-linear decrease of the cotunneling step as a function of mag- 
netic field. This cotunneling step and the non linear evolution could be interpreted as a 
signature of anisotropy [1-5] in N@CgQ . This interpretation relies on the assumption of 
significantly enhanced spin-orbit coupling on the nitrogen atom, which might be caused by 
the hybridization with the gold leads, as the anisotropy of N@Ceo should not exceed tens of 
niT [6-8]. 

The effect of anisotropy can be included in our model adding = —K(S Z ) 2 in our 
Hamiltonian (Eq. 1 of the Letter), where E an i is the anisotropic energy, K is the uniaxial 

5 



anisotropy constant, and S z is the z component of the total spin. However, by including this 
term into the Hamiltonian, we should also observe extra cotunneling steps in the nc 60 = 1 
region. We did not measured this signature, either because of a lack of sensitivity of our 
electronic, or because the finite anisotropy might be charge dependent, but up to now, we 
do not have a clear explanation of the possible charging-state dependence of K. 
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